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We explore consequences of W^-infinity symmetry in the fermionic field theory 
of the c = 1 matrix model. We derive exact Ward identities relating correlation 
functions of the bilocal operator. These identities can be expressed as equations 
satisfied by the effective action of a three dimensional theory and contain non- 
perturbative information about the model. We use these identities to calculate the 
two point function of the bilocal operator in the double scaling limit. We extract 
the operator whose two point correlator has a single pole at an (imaginary) integer 
value of the energy. We then rewrite the Vy-infinity charges in terms of operators in 
the matrix model and use this derive constraints satisfied by the partition function 
of the matrix model with a general time dependent potential. 
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1. Introduction. 



In recent papers [1,2] we have developed a formulation of the c = 1 matrix 
model based on a gauge theory of fermions interacting with a background gauge 
field. This formulation was used to understand the symmetries of the theory as 
well as the general Ward identities relating correlation functions. The natural 
variables in the fermion theory are gauge covariant biolocal operators, y,t) = 
ip(x,t)ip^(y,t), involving fermions at separated points in eigenvalue space. The 
gauge theory formulation led to a natural phase space bosonization of the model in 
terms of variables belonging to a coset space, W± +O0 /H, where H denotes the group 
of transformations which commute with the matrix of eigenvalues of y,t). 

In this paper we explore in greater detail the Ward identities and their conse- 
quences. The main result is an iteration formula for the n-point correlation function 
of the operator W(p, q, t) = 1/2 / dx e ipx ^(x + q/2, t) i/j(x - q/2, t) which deter- 
mines it in terms of all the lower order correlation function. The one point function 
cannot be determined in this fashion and needs to be calculated from the under- 
lying theory, as an input to this system of equations. As a consequence we also 
present an exact equation for the effective action of this theory in terms of a three 
dimensional field. Our equations are completely non-perturbative. As an explicit 
example we calculate the two point correlation function of the bilocal operator ex- 
actly in the double scaling limit. From this general result we also extract the two 
point correlator of the operators W rs = J dx ip*(x, t){(x + id x ) r , (x — id x ) s }tp and 
show that this two point function has a single pole at energy E — i(r — s). Poles at 
discrete imaginary values of energy have been observed in the two point function of 
matrix model operators [3,4,5] and interpreted as signatures of the discrete states 
of the continuum c = 1 theory coupled to gravity. However, the correlation func- 
tions calculated in these works come as a sum of poles, so that the corresponding 
operators couple to a linear combination of these discrete states. Our result thus 
offers a way of isolating the specific operator which couples to a single such state. 

We then proceed to show how the general bilocal operator is related to opera- 
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tors in the original matrix model by using the operator equations of motion. This 
allows us to write down the Ward identities as differential equations satisfied by the 
partition function considered as a functional of a general matrix model potential, 
i.e. Schwinger -Dyson equations of the matrix model itself. We expect that the cor- 
responding differential operators will satisfy a W-jtj.^ algebra. While we have not 
proved this in general, we have checked by explicit computation that this is indeed 
true for the first four equations. These are the "string equations" for this model 
and may be regarded as generalizations of the W p constraints of the p-matrix mod- 
els. These string equations should enable us to understand renormalization group 
flows between various multicritical points. 

2. Gauge Theory formulation of the c = 1 matrix model. 

The double scaled fermion field theory is described by the action 



describes the noncritical string moving in two dimensions Let us introduce the 
operator \^(t) >= J2n c n\ n >> where Cfi clXG operators and {\n >} denotes a 
complete basis in the single particle Hilbert space. The fermion field ijj(x,t) is 
identified with the components of this ket in the x-basis, <x\^(t) >. Introducing 
the background gauge field A(t) with the components A(x, y,t) =< x\A(t)\y > we 
can write (2.1) in the compact notation 





(2.2) 



This action has the background gauge symmetry 



|tf(f) >->W(t)|*(t) > 

A -> U(t)AU\t) + iU(t)d t u\t) 



(2.3) 



where U (t) is a unitary operator in the single particle Hilbert space. 
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We now introduce the operator W(p, q, t) in a "classical" phase space (p, q) by 
the formula 

W(p, q, t) = 1/2 J dx e ipx ^(x + q/2, t) ip(x - q/2, t) (2.4) 

This operator may be also regarded as a generating function for the generators 
\y( r > s ) f the W/^oo symmetry algebra of the critical c = 1 problem [1]. The W^ r,s ^ 
can be obtained from W(p, q, t) by the formula 

= J dp dq F rs (^, P ~jf) W ^ 9. *) ( 2 -5) 

where 

F rs (a,b) = 2cos(ia6) (id a ) r (idt) 8 5® (a, b) (2.6) 

For example 

W^ = (d 2 q -d p ) 2 W(p,q,t)\ p=q = (2.7) 
The operators W(p, q, t) satisfy a simple algebra 

[W(p, q, t), W(p', q', t)\ = i sm(^(pq' - qp')) W{p + p',q + q' , t) (2.8) 

In (2.8) we have explicitly exhibited the string coupling constant g s t r ~ h, so 
that in the limit % — > we regain the structure constants of the algebra of area- 
preserving diffeomorphisms on the plane. From (2.8) we can deduce the algebra of 
the operators W^ r,s ^ defined in (2.5). The symmetries generated by these operators 
have been previously discussed in [1] and [6] in the matrix model and by [7] and [8] 
in the continuum formulation. The quantum deformation of the classical algebra 
of area-preserving diffeomorphisms has been discussed in [9]. Henceforth we shall 
work with h set to unity unless necessary for clarity. 
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One can derive the equation of motion for the operator W(p,q,t) using the 
equation of motion for the fermion field. This equation is 

D t W(p, q, t) = d t W(p, q, t)+ J dp'dq' sm^(pq'-qp') A(p', q', t) W(p+p' : q+q', t) = o| 

(2.9) 

where we have defined 

A(p, q,t) = - J dx e~ ipx A{x + q/2, x - q/2, t) (2. 10) 

The second term in (2.9) can be written as 

J dp'dq' J dp"dq" C^£ l ql A{p\ q', t) W(p", q" , t) (2.11) 

where 

= - P'<l) 6 (P + P'~ P")^ + <t- 9") ( 2 - 12 ) 

is the structure constant of the algebra (2.8). Thus is really the W^-infinity 
covariant derivative with A(p, q, t) as the (background) gauge field. For the critical 
c = 1 theory 

A(p,q,t) = (d 2 q -d 2 p )5(p)5(q) (2.13) 

In this case (2.9) becomes 

[dt+pd q + qd p ]W(p,q,t) = (2.14) 



Equations (2.4) and (2.8) can be understood in a framework that formulates 
quantum mechanics in terms of a classical phase space which is due to Wigner [10]. 
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To explain this, let us introduce the Heisenberg-Weyl group with elements 



9 P , q = e<* p -*V (2.15) 

where (p, q) is a point in the plane R 2 and P and X satisfy the commutation 
relation [P, X] = —ih. The group multiplication law is given by 

9p«9M = e iHpq '- qP ' )/2 9 P+ p', q+q ' (2.16) 
Now noting that g M = e~ iqp '/ 2 e ip * e~ iqp 'I 2 , we see that 

W{p,q,t) = l/2 J dx^(x,t)g M ij(x,t) (2.17) 

One can now immediately see that the algebra (2.8) is an immediate consequence 
of the group law (2.16). 



3. Symmetries 

Among the background gauge transformations of the gauge field A (2.3) there 
are special ones which leave it invariant. These are symmetries of the action (2.2) 
generated by W^ r,s \ Here we ould like to discuss these in terms of the operators 
W(p,q,t) introduced in (2.4). Let us consider the transformations 

6if>(x, t) = --i I dp dq e(p, q, t) e ip ^ x ~^ip(x - q, t) 

If 1 (3 ' 1} 

5ip\x,t) — -i / dp dq e(p,q,t) e tp ^ x+ ^ ip\x + q, t) 

generated on the fermion field by W(p,q,t). Equations (3.1) express the action of 
the Heisenberg-Weyl group on the fermion field. Under these transformations the 



7 



action (2.2) changes by 

SS = j d U M , q M< M A (3.2) 

where 

Dje(p, g, = ^e(p, g, t) + ^ dp'dq' S m^(pq'-qp') A*(p', q', t) e(p+p', q+q', t) (3.3) 

and we have used the hermiticity of A(x,y,t), i.e. A*(x,y,t) = A(y,x,t). Thus 
the symmetries of the action are the set of e(p, q, i)'s which satisfy 

Dje(p,g,t) = (3.4) 

1 theory A*(p,q,t) = A(—p,—q,t) is given by (2.13). In this 

[d t + pd q + qd p }e(p, q,t) = (3.5) 
The corresponding charges are 

W[e] = J dpdq e (p,q,t)W(p,q,t) (3.6) 

where eo(p, t) is a solution of (3.4). It is straightforward to check that (3.4) is sat- 
isfied by e(p,q,t) = e" (r-s) *F ra (2^, thus verifying that the generate 
the symmetry algebra in the critical c = 1 model. 



For the critical c = 
case (3.4) reads 
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4. Ward Identities and the three dimensional effective action. 

For more general e(p,q,t), not necessarily satisfying (3.4), the variation of the 
action in (3.2) leads to Ward identities for the correlation functions of W(p,q,t). 
We now turn to a discussion of these. 

Consider a general n-point function of the W(p,q,tys 

G n (Pi, ft, U) =< /i\T(W(p h q h h), W(p 2 , q 2 , h) • • • W(p n , q n , t n ))\n > (4.1) 

where by \fi > we have denoted the filled fermi vacuum. In the operator formalism 
the Ward identities can be derived by using the equation of motion for W(p, q, t), 
(2.9), together with the symmetry algebra (2.8). For example, for the two point 
function we have 

D t G 2 (p,q,t;p',q',t') = D t 0(t - t') < »\W(p, q,t)W(jpf , q' > + 

D t 9(t' -t)< »\W(j/, q', t')W(p, q, t) |/i > 

Using the equation of motion (2.9) and the symmetry algebra (setting h — 1) we 
get 

D t G 2 (p, q, t; p', q', t') = S(t - t') < fx\ [W(p, q, t), W(p', q' , t')\ \» > 

1 (4-2) 
= i6(t - t')sm-(pq' -p'q) Gi(p + p , q + q ,t) 

Alternatively one may use the functional integral formalism to derive the Ward 
identities. In the functional formalism the n-point function (4.1) is given by the 
expression 

If n 
G n (pi, q h U) = z V^V^e tS J] W(pi, q u U) (4.3) 

where Z is the partition function. To derive the Ward identites we make the change 
of variables ip — > ip + dip, ^ — > ift + 5ip^ where the dip and 8ip^ are given in (3.1). 
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Under this change of variables, the action S changes as in (3.2) while the change 
in W(p, q, t) is 

5W(p, q,t) = j dp'dq' siv^{pq' - qp') e(p' , q' , t')W(p + p',q + </, t) (4.4) 
Thus we get 

1 r n 
= - / ThfiThl> e lS [i5S Y[W(p t , qi ,U) + 

A J i=i 

n 

^2 W(Pl,?l,*l) • ••8W(pi,q i ,ti) ■ ••W{p n ,q n ,t n )] 
i=l 

(4.5) 

Substituting (3.2) and (4.4) in (4.5) we get the Ward identity 
B t G n+ i(p,q,t;pi,qi,ti) 

n 1 

= i^2S(t- tj) sin-ipqj - qpj) G n (p h qi, h; ■ ■ ■ ; Pj + p, qj + q, tj; ■■■;p n , q n , t n ) 

(4.6) 

Equation (4.6) gives a closed set of first order differential equations for the 
multipoint correlation function of the W(p, q, i)'s. It is appropriate to point out 
that this has happened because not only is the change in W(p,q,t) under (3.1) 
proportional to another W operator, but more importantly the variation of the 
action also involves a W(p,q,t) (equation (3.2)). The situation may be compared 
with that for the correlation functions of the Wilson loop operators in large- iV QCD 
where a closed set of Dyson-Schwinger equations can be derived precisely because 
the variation of both the Wilson loop operator as well as the action involves another 
Wilson loop operator. Starting with the one point function which can be calculated 
in the fermion field theory (2.1), equation (4.6) can be used repeatedly to solve for 
the n-point function. 

These Ward identities can be summarized in terms of a differential equation 
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satisfied by the generating functional 



Z[J(p,q,t);A] = J V^Vtpexp [iS + J dpdqdt J(p,q,t)W(p,q,t)] (4.7) 

The equation is 

S f 1 & 

[^777 — + / d p' d( i' sin o(p'i - <i'p) j (p - P, ?' - 9, r 7 , , , tt ]^; A] = 

oJ(p,q,t) J 2 5J(p',q',t) 

(4.8) 

The above equations lead to an equation satisfied by the effective action defined 
by a Legendre transformation 

r[$(p,q,t);A] = -log Z[J; A] + J dp dq dt J(p,q,t)<f>(p,q,t) (4.9) 



Since 
we have 

/l (5 
dp'dg' sm-(p'g - q'p) $(p + p', g + g', t) ^ - - r[$; A] = 

(4.11) 

Equation (4.11) gives an exact formulation of the c = 1 matrix model in terms of 
a three dimensional field q,t) . 

For the critical c = 1 model (4.6) reduces to 



[dt+pd q + gcyG„+i(pi,gi,t;pi,gj,tj) 

n 

% 



^2S(t-tj) sin-(pqj -qpj) G n (pi, qi, h; ■ ■ ■ ;pj +p, qj + q, ty, ■ ■ ■ ;p n , q n , t n ) 
3=1 

(4.12) 

Because of time translation invariance, in this case the one point function G\ is 



* A three dimensional action for the discrete states of the c = 1 Liouvillc theory has been 
proposed in [8]. 
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independent of time, so that we have the equation for G\ 

(pd q + qd p )G 1 (p,q,t) = Q (4.13) 

This implies that G\(p,q,t) is a function of (p 2 — q 2 ) only, as may be explicitly 
verified by a computation of G\ in the critical fermion field theory. 

Before ending this general discussion on Ward identities we would like to men- 
tion that we have so far obtained identities satisfied by correlation functions of the 
generators of VT-infinity transformations of the theory. There are other operators 
in the fermion field theory besides these, viz. operators which transform homo- 
geneously under W^-infinity transformations. Examples of such operators are the 
group elements of ^-infinity. There would be Ward identities satisfied by corre- 
lation functions of such operators as well. The situation is analogous to WZNW 
models where apart from the currents there are the primary fields g{z,z) which 
are group elements and there are the Knizhnik-Zamolodchikov equations involving 
the correlators of these primary fields^. 



5. The two-point function 

In this section we focus our attention on the Ward identities for the c = 1 
critical theory. Of particular interest is the two point function G2(p,q,t;p',q',t') 
of the W(p, q, t)'s since this contains information about the spectrum of the theory. 
In this case, equation (4.6) gives 

[d t +pd q + qd p }G 2 (p,q,t;p',q',t') =i5{t-t ! ) sm^(pq'-qp') G^p + p' ,q + q') (5.1) 

where we have used the fact that G\ is independent of time. In terms of the fourier 
f We would like to thank S. Shatashvili and E. Verlinde for a discussion of this point. 
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transforms 

W(p,q,E) = J dt e~ iEt W(p,q,t) (5.2) 
equation (5.1) can be written as 

[iE + pd q + qd p ]G 2 (E;p,q;p', q') = 2m sin-(pq' - qp') d(p + p' ,q + q') (5.3) 

where we have defined G 2 by 

J dt dt 1 e~ iEt e~ iE '* G 2 (p, q, t; p'q', t') = 5(E + E') G 2 (E; p, q, ; p', q') 

Note that since G 2 is symmetric in the interchange of the arguments (p, q, t) — > 
(p',q',t'), G 2 has the symmetry G 2 (E;p, q, ;p' , q') = G 2 (-E;p'q';p,q). 

The zero modes of the operator [iE + pd q + qd p ] are precisely the solutions 
of the equations of motion. They occur at E rs = i(r — s) and are given by the 
functions ^rs(^^,^^) which were introduced in (2.6). 

In order to solve (5.3) we continue the energy E to complex values. Absorbing 
the i in the definition of E we have 

G 2 (E;p,q;p',q') = E + J qQ [2ttz sm±(pq'-qp') G^p+p'f -(q+q'f)} (5-4) 

Using the exponential representation of the inverse operator we get 

(sgn)xoo 

G 2 (E;p,q,;p',q') = J ds e~ sE e ~ s ^ + ^ [2m (sin|) G 1 (R + R' + 2P)} 



(5.5) 

where we have defined 

P = pp' - qq' Q=pq' - p'q R = p 2 - q 2 R' = p' 2 - q' 2 (5.6) 
Now, noting that the operator (pd q + qd p ) is the generator of Lorentz rotations in 
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the (p, q) space we can easily write (5.5) as 

(sgnE)xoo 

G2(E; p, q; p', q) = 2?ri J ds e~ sE sm[-(Q cosh s — P sink s)] 



x Gi(R + R' + 2P cosh s - 2Q sink s) 

(5.7) 

The solution (5.7) satisfies the symmetry requirement G(E; p, q, ; p', q') = G(—E; p', q' 
This expression for G2 exactly matches with the explicit calculation of the two point 
function for q = q' = given in [4]. 

In a similar way one can obtain the three point function in terms of the two 
point function and the one point function and so on. 

To proceed further we need to know the one point function G\. This can be 
computed in the fermionic theory. Using the notations and conventions of [4] * the 
result may be written as 

00 

Gi(p,q) = ^ J dxe ipx J dv ^ VeO, x + ^q)M", x - ^q) (5.8) 

Taking a derivative of (5.8) with respect to \i we get 

1 f \ \ 

d^Gi{p, q) = -- J dx e ipx M»> x + ^)M^, ^ ' gS) ( 5 - 9 ) 



e=± 

Using the relation (A. 12) in [4] one can show that 

o^+i(p 2 -9 2 )coth| 

4 



d^G^q) = -^Re / di — (5.10) 

4v 2tt J smh§ 



The integral on the right hand side is well defined for complex (p 2 — q 2 ) with a 
small positive imaginary part. The result is then analytically continued for real 



* Note that our x stands for the matrix eigenvalue while time is denoted by t. Also, the 
eigenvalue variable A used in [4] is y/2x. 
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(p 2 — q 2 ). Now, taking a /i-derivative of (5.7) and using (5.10) one gets 

oo (sgn E)xoo 

d fJ/ G2(E;p, q;p', q') — =Re / j / ds e" s£; sin-(— Psinh s + Qcosh s) 

2y2 J sinh| J 2 



2 



x exp [ifj£ + -(R + R' + 2Pcosh s - 2Qsinh s)coth|] 

(5.11) 



or, 



OO 

I r e i^+i(R+R')coth 

d /Jl G2(E;p, q, ;p', q') = -— =Im / d£ 



4^2 7 sinhf 

o z 

{sgaE) Xoo 



, . ^Pcosh(s — |) iQsinh(s — |) 



x / dse- sl "[exp( y — r ^-) 

J 2sinhi 2sinh| 

ll 

?Pcosh(s + I) iQsinh(s + §) 

exp ( t t )] 

2sinh| 2sinh| 

(5.12) 

For q = q' = this expression for d^G^ matches exactly with that given in [4]. As 
in that reference, the integration over s may be performed by changing variables 
from s to y = e s and expanding in inverse powers of y. The result is most naturally 
expressed in terms of the new variables z±, z'± defined as 

z± = ^=(p±q) 4 = -^(p'±</) ( 5 - 13 ) 

In terms of these variables we finally get 

oo 

d^G 2 (E; p, q, ; p', q') = 3-=Im f e ^+|(* + *-+4*-)coth§ 

o 2 



,Z+z'_ , E 



- ' — i — - . tt|bi sinh \E\i , d 'z+Z-zLz'_ . 

2-2+ sin tt|P| 1 1 S i n h § 



00 / i / i <>• / 7 7~~ 

E nr/ Z+Z_.n 1 n 1 - , C . \/ Z-L Z— Z i Z_ 
i n {(—^)~ t =r + (-7^)» 7 ftt} sinh JJ^—^ 7—)} 
x n z_z' + J (n-E) K zLz+ J (n + E) S 2 n{ sinh i n 

(5.14) 
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Not surprisingly, equation (5.14) shows more structure than the corresponding 
expression for q — q' — 0. In particular, we can now obtain the two point func- 
tions of the operators W rs (E), defined in terms of W(p, q, E) (eqn. (5.2)) by an 
expression similar to (2.5) 

W rs (E) = 2(-id^n-id + yicos(z_z + /2)W(p,q,E)} z±=0 (5.15) 

where d± = The limit z± — > in (5.15) has to be taken carefully. In fact, a 
sensible limit is defined by going to a region of the complex E plane in which the 
expressions are well-defined. To see how this comes about, let us compute the two 
point function < W rs (E)W(p', q', — E) >. We will assume that Re E is positive 
and that r < s. The other cases can be treated similarly. This two point function 
can be projected out from (5.14) by taking suitable derivatives with respect to z± 
as given in the definition of W rs (E) in (5.15). We will calculate only the leading 
behavior for large fx, which corresponds to the planar limit. In this limit the terms 
in which the derivatives in (5.15) hit the cosine factor are lower order in fx. So we 
can set this factor to unity and directly take the derivatives d r _ d b + on G2 in (5.14). 
To perform the limit z± — > we first work in the region Re E > (s — r) where 
the limit exists and then analytically continue to other regions of the complex E 
plane . 

Let us now compute the two-point function < W rs (E) W r > s >(—E) >. As before 
we will assume that E is positive and consider the case r < s. It turns out that 
the above two point function vanishes unless (s — r) = (r' — s'). Thus the region 
in complex E plane where the limit z± — > has to be performed in order to define 
W rs (E) is precisely the region where W r / S i(—E) can be defined. We thus have a 
completely well defined two point function. In fact, we find the leading fx behavior 



* The limit E — > Z is more tricky and one needs to be careful while taking this limit. In fact, 
just as for q = q' = 0, the poles in the expression in square brackets in (5.14) cancel in this 
limit. 
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to be 



< W rs {E)W r , s ,{-E) >~ 5 s - r y^^ r+s+r ' +s ^ log n r [S . T) (5.16) 

h, - [s - r) 

It is interesting to note that the above two point function contains a pole 
at a single integer value of the energy, which is actually a consequence of the 
fact that W rs are eigenoperators of the hamiltonian. In fact a standard matrix 
model operator O m = f dx x™ 1 ^^ ma y De written as binomial sum of W rs 's with 
r + s = m. This explains why the two point function < O m O n > contains 
a sum over poles : the operators O m are not energy eigenoperators, but linear 
combinations of energy eigenoperators. 

It is also interesting to note that the above expression for the two point function 
is exactly reproduced in the analytic continuation of the theory described in [1]. 
In the phase space formulation one continues t — > it as well as p — > — ip (where p 
denotes the single particle momentum). Note that this analytic continuation is not 
the standard euclidean continuation. It leads to a hamiltonian h — —\{p 2 + x 2 ). 
The fermi sea is filled from — oo to some level —fi. The operators W rs are well- 
defined in this framework and their two point function may be computed in a 
standard manner. In the planar limit the result is exactly (5.16). However, in 
this framework the "tachyon" pole corresponding to the first term in (5.14) is not 
obtained since the energy levels are all discrete. The equivalence of these two cal- 
culations is not entirely unexpected since the calculation of two point functions of 
operators like J dx x m ^ip using harmonic oscillator wavefunctions with imaginary 
frequencies [5] agrees with that of [4]. 

Finally we would like to emphasize that in the Minkowski section, E in (5.16) 
is purely imaginary so that in Minkowski space the poles are at purely imaginary 
values of the physical energy. Consequently they cannot be associated with states 
in the spectrum of the theory. 
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6. c = 1 Matrix Model with a general 
potential and String equations 

At first sight the operators W(p,q,t) for q ^ do not have an obvious corre- 
spondence with operators of the original matrix model. This is because the most 
general invariant local operator in the matrix model is of the form Tr exp(ikM(t)) 
whose expression in the fermion field theory is simply W(k,0,t). In this section 
we will show that the q ^ operators can be directly expressed as matrix model 
operators through the Heisenberg equations of motion. 

Consider a matrix model with an arbitrary time dependent potential U(M(t)). 
The corresponding fermion field theory (2.1) has a background gauge field 

A(x,y,t) = ^(d 2 x -V(x,t)) (6.1) 

where V(x, t) denotes the corresponding scaled potential. Corresponding to this 
gauge field the expression for A(p,q,t) defined in (2.10) is 

A(p, q, t) = (d 2 q - V(id p , t))5(p) 5(q) (6.2) 

Therefore, in this case the equation of motion (2.9) for W(p,q,t) reads 

[9t + Pd q - \[V(-id p + lg, f) - V(-id p - ±q, t)]]W(p, q,t) = (6.3) 

This equation may be used to relate the operators 

W n (p,t) = d2W(p,q,t)\ q=0 (6.4) 

to matrix model operators O n (t) = TrM n (t) = J dx x n ^ (x,t)ip(x,t). Using the 
above definition of W n (p,t) we can obtain, from (6.3) 

W,0+P^+i(p,t)-it|^Q V<^ m X-id p ,t)W m (p,t) = 0, n = 0,l, 

m=0 

(6.5) 

where = 1 for k odd and zero otherwise and we have used the notation V^ k \x, t) 
to denote d^V(x,t). Clearly we can use (6.5) repeatedly to relate W n (p,t) to 
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Wo(p, t) by an equation which contains at most n derivatives with respect to time. 
The Wo(p,t) may be used, in turn, to generate all the operators O n (t). In fact, 



O n (t) = 2(-id p ) n W (p,t)\ p=0 (6.6) 

(The factor 2 comes from the definition of W(p, q, t which has a factor of \ in it). 

The above relations mean that the set of operators W n (jp,t) which have no 
direct correspondence with the operators of the matrix model, is equivalent to the 
set of operators d™Wo(p,t) which are simply <9«-Tre ipM ^ . The role of q is replaced 
by time derivatives. 

It is now clear that the equations (6.5) imply an infinite set of equations for the 
one point functions of the O n (t); the equation for O n (t) involving at most (n + 1) 
time derivatives. Equivalently, these equations may be expressed as an infinite set 
of constraints on the partition function, which we may write as 

W n (t)Z = 0, n = 0,l,2,-.. (6.7) 

These may be written in terms of the time dependent couplings V^(t) defined as 

oo 

V(x,t) = Y,V k (t)x k (6.8) 
fc=i 

Below we list the first few of these constraints 

k — 1 

mt) = a > wt) + + mm WD + kd ' VM mr> ] (eAl) 

k — 1 



19 



W3(i) " Sf + 2 y (t + i) nm iwi + 'ivi+iW 

3, ,2,.,., * 3 * 



q 00 A 



4^ -^WiW 

(6.12) 

Alternatively, these equations may be derived by performing a set of specific 
W transformations with time dependent coefficients on the fermionic fields, in such 
a way that the transformed action is of the same form as the original action but 
with a modified potential. These transformations are generated by the charges 

n „ 

I n = J2 dx ^(x,t)e k (x,t)[-id x } k ij(x,t) (6.13) 
where the functions e^it) satisfy the following set of equations 



j -y n-m , . ^ ^\ 

id t €m + 2 d x e ™ - g € j+m(-iy +1 i J y ^ Jd J x V(x,t) + i^Cm-i = 



(6.14) 



m — 1, 2, • • • (n — 1) 



da;e n _i + <%e n = 0, e n = e n (t) 

Clearly these equations can be used to determine any ej in terms of e J+ i and thus 
express all of them in terms of e n . Since e n (t) is a function of t alone we can choose 
a fourier basis to write the constraints on on the partition function as 



6 Z = 01 



W n {E)Z =[j 'dxdte- i ^[^eo+^e -^e J -(x,t)(-iy + ¥^V(x,t)] (yr( ; 

(6.15) 

The operators W clearly form a closed algebra. Since the equations (6.15) are 
rewritings of the original Ward identities, this algebra is expected to be isomorphic 
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to the W - algebra . The role of spin is played by n, while the "mode number" 
is E. Note that the mode number can be any real number and not necessarily an 
integer. 

For time independent potential and for time independent e n in (6.13), I n re- 
duce to the well known transformations which generate the isospectral flows of the 
Schrodinger operator 

lk 0) = (-d 2 x + V(x))%\ + (6.16) 

where |+ denotes the part of (6.16) which involves only positive powers of the 
differential operator — id x . and the corresponding string equations are simply given 
by 

/ dxd x R n (V(x))-J--Z = (6.17) 
J 6V(x) 

where R n (V(x)) are the Gelfand-Dickii coefficients [12]. For the special case of 
time independent quadratic critical potential V(x) = —x 2 these isospectral defor- 
mations have been discussed in [11]. Time dependent transformations I n and the 
corresponding algebra for the special case of V(x) = —x 2 have been given in [6]. 
A different set of identities (once again for time independent potentials) on the 
puncture operator, rather than the partition function, have been derived in [5]. 

The W n equations would be probably useful in studying the flows between 
various multicritical points. In understanding such flows one must remember that 
the potential V(x, t) used above is a scaled potential. In general, starting from some 
matrix model potential U(M) in which the N dependence of the various terms 
appear in such a way that the large- N expansion is a topological expansion, the 
critical potential in the double scaling limit will contain only a single monomial in x. 
To obtain a scaled potential V(x) which has various powers of x the N dependence 
of the terms in the original matrix model potential have to be different from the 
standard N dependence. This does not, however, prevent us from discussing flows 



★ We have explicitly checked this for n = 0, 1, 2, 3 
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between various multicritical points. Rather, this means that when we consider, 
e.g. flows between a critical point of the form V(x) = —x 2 and another critical 
point V(x) = —x 4 by an interpolating potential V(x) = ax 2 + 6x 4 , only the end 
points a = —l,b = and a = 0,6 = — 1 correspond to random surface theories. 
For generic (a, b) the matrix model jj expansion no longer correspond to a string 
perturbation theory. 



Note added : After this work was complete we learnt that I. Klebanov and A. 
M. Polyakov [13], D. Kutasov, E. Matrinec and N. Seiberg [14] and G. Moore, M. 
Plesser and S. Ramgoolam [15] have also considered the problem of constraining 
correlation functions using different methods. 
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